Generalizations of Kijowski's time-of-arrival distribution for interaction potentials 
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Several proposals for a time-of-arrival distribution of ensembles of independent quantum particles 
subject to an external interaction potential are compared making use of the "crossing state" concept. 
It is shown that only one of them has the properties expected for a classical distribution in the 
classical limit. The comparison is illustrated numerically with a collision of a Gaussian wave packet 
with an opaque square barrier. 
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. I. INTRODUCTION 

X) . 

^ ' In many experiments, the observed quantities are the instants of occurrence of certain events, or the durations of 
processes. However, the standard quantum formahsm does not provide, at least in an obvious manner, a working rule 
for time observables. In spite of the difhcuhies and objections to consider time as a quantum observable (by Pauli 
[01, AUcock Q, and other authors), many researchers have sporadically tried to fill this theoretical lacuna. The effort 
in that direction has become more intense and systematic in recent years. In particular, much attention has been 
devoted to the quantum description of the "arrival time" , see for a recent review. In its simplest form, the 

problem is to define an ideal quantum arrival time distribution for a wave packet moving freely in one dimension. By 
imposing a number of classically motivated conditions (normalization, positivity, minimum variance, and a certain 
CnI , symmetry with respect to the arrival point X) this is solved uniquely by Kijowski's distribution [^. This distribution 
■ may also be associated with the positive operator valued measure (POVM) generated by the improper eigenstates of 
' the time-of-arrival maximally symmetric operator of Aharonov and Bohm p^]. Tab, 
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' where x and p are position and momentum operators and m is the mass. Aside from the discussion of remaining 
J2 > interpretational puzzles, an important pending question is its generalization for particles affected by interaction 
potentials However, Kijowski's set of conditions cannot be applied in the general case ||^, where, classically, not 
all particles necessarily arrive, or multiple arrivals may have to be considered. Neither can the simple symmetrization 
rule leading to Tab be used pof . Thus a generalization of Wk required some novel approach and it is only very recently 
that some of them have been explored, following different heuristic and/or formal arguments. Their physical content 
; ^ must be analyzed in order to select one that does fit into the proposed objective (although several might be adequate). 
I This is the aim of the present paper, where we shall examine three possible generalizations of Kijowski's distribution. 



The unifying framework is provided by the concept of "crossing state" introduced in |23|, and inspired by Wigner's 
formalization of the time-energy uncertainty principle [ pO| . In all the proposals examined here the (candidate) time- 
of-arrival distribution at point X may be obtained from the overlap of the time dependent wave function and the 
crossing states \u^^{X)) as 

H(T) ^Y.^P^T) = mT)\u^X))\\ (2) 

where /? = L, i? is an index for "left" and "right" ; its exact meaning will vary in the different generalizations. We shall 
use a unified notation that will differ in general from that of the original papers to facilitate the common presentation 
and comparison. Incidentally, all distributions defined in this manner are automatically covariant with respect to 
time translations, namely, the arrivals predicted for a given fixed instant are independent of the choice made for the 
origin of time p5|,pl. 



II. DIFFERENT CROSSING STATES 
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A. Free motion 



Kijowski's distribution is obtained from the general expression (|]) with the foUowing crossing states (their coordinate 
representation, time evolution, and wave packets peaked around them have been studied in 

\u^{X))k = {\p\lmf/^Q{ap)\X), (3) 

where 

- - if/3 = L 

- if /3 = i?, ^ > 

\p\^^'^ is defined by its action on plane waves, 

\p\^/^\p)^\p\^/^\p), (5) 

(the positive root is taken), and Q is the Hcavisidc distribution. 
Now we can write Kijowski's distribution in operator form as 

^k{T) = [0(«P)(|p|/m)i/2<5(X - xm/mf'Mo^p)] |^(T)). (6) 

a 

Each of the operators in brackets (with a = + or — respectively) corresponds classically, i.e., disregarding the lack of 
commutativity between x and p, to the classical dynamical variable 

5{x~X)^Q{ap), (7) 
m 

whose average {■■■■) over a classical phase space density gives the flux due to particles arriving from the left, J^, or 
minus the flux due to particles arriving from the right, (which is a positive quantity, > 0), 



a 



jj^^^{S{x~X)-^Q{ap)),„ (8) 



in other words, the modulus of the flux of particles of the classical ensemble that arrive from one side at a given time 
. The addition of these two contributions is the classical time-of-arrival distribution, 

nc/(r) = - jS (9) 



B. Interacting case: first proposal 

The previous discussion motivates the first generalization of Hk considered here for independent particles affected 
by an arbitrary interaction potential. Since (|9|) is valid regardless of the presence of a potential, and the dynamical 
variables for the two flux contributions are always given by (0), it was proposed in that the quantum time-of- 
arrival distribution in the general case be given by the same expression used in the free motion case, Eq. (|^), and by 
the same crossing states, \u^)i = \u^)j^, as they lead to operators in correspondence with the classical expression (|^). 
In other words, the proposed distribution is given by 

ni(r) = ^|(^(r)|«%p, (lo) 

where now ip evolves with the full Hamiltonian H = Hq + V which includes a kinetic term, Hq, and a potential term, 
V. We are thus extending to the quantum domain the fact that in classical mechanics the expressions for dynamical 
variables representing the partial fluxes do not vary from the free-motion case to the interacting case. Here, the 
state (ipiT) quantally, and the evolved phase space density classically), rather than the dynamical variable, contains 
the information which is specific to each particular Hamiltonian. At the very least, this generalization of Kijowski's 
distribution has the merit of being simple, arguably the simplest one. Further properties were commented in the 
original paper psf . In particular, note that ( |lO[ ) need not be normalized, and in fact may be not normalizable, as may 
also be the case classically, e.g. because of periodic crossings in a harmonic potential. It can also be defined even if 
the system is not classically integrable. 
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C. Interacting case: second proposal 



The other two definitions discussed here are apphcable to "scattering potentials" where MoUer operators exist, 

n± = hm e^Ht/h^-^Hot/h^ (^^^ 

t — >=poo 

We shah assume that the potential operator V has a local coordinate representation 

{x\V\x') ^ 6{x ~ x')V{x), (12) 

with potential function V{x) vanishing at long distances from the interaction region so that the (strong) limits in Eq. 
( pi] ) exist. As usual, we shall also consider, with the same notation, extensions of these Hilbert space operators that 
can be applied to plane waves, 

\P±) = f^±b) = \P) + Im ^V\p), Ep = pV2m. (13) 

e-^o+ Ep±is ~ H 

The Lippmann-Schwinger states \p±) are improper (not square integrable) eigenstates of H with eigenvalue Ep. fij- 
are in general only isometric, i.e., they conserve norm, but in the absence of bound states they become also unitary. 
We shall limit ourselves to this later case hereafter. The physical meaning of n± is best understood with the aid 
of "asymptotic" incoming and outgoing states, and (f>out, to which the actual state ip tends in a strong sense in 
the infinite past and future, respectively. These asymptotic states evolve freely, with Hq. Whereas 17+ provides the 
scattering state by acting on the incoming asymptote, r2_ does the same job acting on the outgoing one, 

f2_|</)o„t(t))-|^(t)}, (15) 

for all t. 

Without bound states, the resolution of unity can be written equivalently as 

/oo poo 
dp\p){p\= dp\p±)(p±\, (16) 
-oo J —oo 

and the Hamiltonians Hq and H may be related by the unitary transformation 

H = h±Hah±, (17) 
with a similar relation holding for functions of Hq and H. This suggests using the crossing states 

\ui{X)),^h±\u^'{X))j,. (18) 

Note the additional subscript in the crosssing state, + or — , due to the possibility to act with either one of the two 
MoUer operators on \u^)j^. Thus, this procedure generates two different distributions labelled by + or — , 

Il,^±=Y.\mT)\uiiX)),\\ (19) 

These distributions previously appeared in |po| (later superseded by [pT| ). It is important to notice that this is not 
the way in which the original distributions in ]20[ were presented; they have been adapted here to the unified notation 
we are using in order to achieve a better handle for comparison. Their physical content is made evident by making 
use of (jl^), (p^), and the isometry of the MoUer operators, 

{i^{T)\n+\u^{x))j, - {^UT)\nlh+\u^{x))j, = {^ut)\u^{x))k, (20) 
{^{T)\n_\u0{x))j, = {cbout{TmlhM'{x))K = (0o«t(r)|w^(x)}^. (21) 

This means that the generated distributions, 112^+ and n2^_ are nothing but the Kijowski distributions corresponding 
to the incoming and outgoing free-motion asymptotes, respectively. These may be useful objects before (112.+) and 
after (112, _) the collision, but not in the midst of it. 
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D. Interacting case: Third proposal 



A different, more complex proposal is based on the following crossing states, 

|4(X))3 = ^±Q{amP\I^Y'^^i\X), 
which lead again to two different distributions, 

n3,±(r) = ^|(V'(T)|4(x))3p. 





(22) 



(23) 



These distributions were first introduced and discussed by Leon et al. in |2] 
quantum canonical transformations of the classical time of arrival (see also p2 | 
As previously stated, the rewriting in terms of crossing states is intended to clarify the physical consequences of 
these distributions. To analyze their meaning, let us first study the amplitudes corresponding to ri+, by inserting 



as justified by quantization through 
for further analysis of this proposal). 



1 = fJ^^+j ^-iid rewriting them as 



(VXr)i \n+Q(aml] \?l+{\p\/mf/^?ll] \X). 



(24) 



The operators involved have been separated in two brackets that can be interpreted physically. The first one is a 
projector that selects the part of a wave function that had positive (a — +) or negative momentum (a = — ) in the 
infinite past. 



f1 = n+e{ap)n^ 

the second group of operators is 

(2H/m)i/4 = m-i/2f)+|p|i/2f}t 



t _ 



dp\p+){2E,/mY'\p+\ 



(25) 



(26) 



where the positive root is taken. For the case in which the incoming asymptote is restricted to positive momenta, 

(^(T)|F^ = (V'(T)|, (27) 
(^(r)|F^^-0, (28) 

so that the proposed distribution takes the form 

n3,+ (r) - (T) - {i,{T)\{2H/mY/H{x - X){2H/my^^mT)). (29) 

The resulting operator is a quantum symmctrization of the classical phase space dynamical variable 

{2E/my/^dix-X), (30) 

E being the total energy, to be compared with the classical variables of the first generalization which lead to the 
positive and minus negative fluxes. 



{\p\/m)Q{p)S{x - X), 
{\p\/m)e{-p)S{x-X). 



(31) 
(32) 



Note that, for an X such that V{X) ~ 0, the (classical) average of (|3^) is equal to the average of the sum of ( |3l| ) and 
(^), i.e. to the time of arrival distribution. 



{{2E/m)'/^S{x ~ X))^, = Jii - , if V{X) = 0. 



(33) 



However, in the interaction region, it is in fact proportional to the local square root of the total energy. This may 
lead to significant differences with the actual time-of-arrival distribution, which classically is always given by J^^ — 
irrespective of the value of V{X). In particular, for a potential barrier such that the initial (asymptotic) momenta 
of the particles slow down on its top to a smaller value, (^0|) will overestimate the value of the true arrival time 
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distribution, since for each trajectory the square root {2mEy^'^ is used, instead of the smaller local momentum that 
determines the partial fluxes. Note also that the superscript L in (^9|) does not make reference here to actual crossing 
from the left in the classical limit, but to motion from the left in the infinite past. In this respect, (^) equally counts 
right or left crossings at a given time. 

A similar analysis may be carried out for a state with incoming asymptote in the subspace of negative momenta. 
Thus, for an arbitrary initial state with positive and negative momentum components, 113^+ = ^3 + + 11;^^ provides 
a quantum version of the distribution that corresponds to the classical dynamical variable (^0|). Repeating the same 
steps with Ilg _ may be interpreted as the quantum versions of the classical distributions of ( ^0[) for particles 
that will have positive or negative momenta in the infinite future. 



III. NUMERICAL EXAMPLE 



We shall illustrate the results of the previous section with the collision of a wave packet for a particle of mass 
m — 1 with a square barrier. The initial state (at i = 0) is chosen as a minimum-uncertainty-product Gaussian with 
negligible negative momentum components. Its average position, average momentum, and standard deviation are for 
the first four figures, in atomic units, 

(S) = -6, (p)=6, Ax^l, (34) 

whereas the barrier energy, initial and final points are 

Vn = 10, X, = 0, Xf ^ 10. (35) 

With these parameters the barrier is rather opaque. This prevents tunneling and makes the collision predominantly 
classical. The figures combine variously the different distributions, their components, and the flux for the three points, 
X = —2, 5, 12, corresponding to positions before, in, and after the barrier respectively. The evaluation of the integrals 
requires explicit representations of the states \p±) that can be found elsewhere 

3 I ■ < ' < ■ < ■ < ' 1 




T (a.u.) 

FIG. 1. Hk and 112,+ (which overlap) for X = —2, 5 and 12 (solid, dotted, and dashed lines) for the Gaussian wave packet 
described in the text, see 

Fig. 1 represents Kijowski's distribution for the case in which the wave packet evolves freely, without interaction 
with the barrier, as well as 112.+ in the interacting case. Since the fraction of negative momenta in the initial state 
is negligible, it may be considered essentially equal to the corresponding incoming asymptotic state for the numerical 
accuracy of the flgures. Thus 112.+ (or 112 ^) ^^"^ are indistinguishable for the three values of X. 

Fig. 2 shows IIj _ and the flux J, before, in, and after the barrier, which is now present. They coincide after the 
barrier, for X = 12, but differ otherwise. This is to be expected since 112 _ a free-motion Kijowski distribution for 
the positive momentum part of the outgoing asymptotic state. 
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T(a.u.) 

FIG. 2. and J (solid and dashed lines respectively) for X = —2,5 and 12 (indistinguishable in the later case). 

At X = —2, where the potential vanishes, Hi and 113^+ arc in essential agreement, see Fig. 3. The second, smaller 
bump corresponds to the period of negative flux due to a small reflection. This smaller bump is due to a contribution 
of Ilf only, with no contribution from Ilf , since this negative flux is associated with crossing from the right. However, 
^ alone (which in the figure is indistinguishable from the total Hs ^) provides incident and reflected bumps even 
though they are associated with different crossing directions. Recall in this respect that the superscript L in Hg _|_ 
means "motion from the left (positive momentum) in the infinite past" . (The reflected part will be seen in more detail 
in a different collision described in Fig. 5 below.) The difference nf — Ilf has also been depicted; note its agreement 
with the flux. 

3 I ■ < ' < ■ 1 




-1 I ■ ' ' ' ■ 1 

12 3 

T (a.u) 

FIG. 3. J (dots), n3,+ (solid line). Hi (long dashed line), and nf - nf (short dashed line) for X = -2. 

At X = 12 there is only transmission and Hi, n3_+ and J are indistinguishable. However, for X = 5, i.e. in 
the barrier, H3^+ is clearly larger than J and Hi, which essentially coincide, see Fig. 4. A simple estimate of 

the ratio between the peaks follows from the classical limit discussed in the previous section: fli must be to a good 
approximation proportional to an average local momentum; with the current parameters this entails [2(18—10)] = 4. 
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It is important to notice that the relevant quantity is the total energy (taken as w 18 from the initial value of the 
momentum) minus the potential barrier energy. On the other hand, n3^_|_ is proportional to the square root of the 
total energy, (2 x 18)^/^ = 6. The ratio of the two peaks in the figure is indeed 6:4. 

2 I ' . ' . ' 1 




-0.5 I • ^ ■ ^ • 1 

13 5 7 

T (a.u) 

FIG. 4. J (dots), n3,+ (solid line), and Hi (dashed line) for X = 5. 

As a complement to Figure 3, we have lowered the average momentum of the initial wave packet to (p) = 3, while 
keeping Ax = 1, so that the whole packet is now reflected. We evaluate J and all distributions at a; = — 5 for an 

initial average position (a;) = —9. The flux, as portrayed in Fig. 5, shows clearly a positive part during the incidence 
and a negative part corresponding to reflection. Since several combinantions of the different distributions and their 
components match with adequately the two bumps of |J| or just one of those, we have in fact only represented 
113^+ = + , with a dashed line. IIi is barely distinguishable from it. The incidence bump on the left is also 
reproduced by 112,+, or by Ilf , whereas the reflection bump is reproduced by 112,- or H^. 

1.5 , ■ ^ ■ ^ ■ ^ ■ ^ 1 




-1 I ■ ^ ■ ^ ■ ^ ■ ^ 1 

2 4 6 8 

T (a.u) 

FIG. 5. J (solid line) and 113,+ (dashed line) for X = —5. See the text for details. 
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IV. CONCLUSION 



We have examined three different generalizations for Kijowski's time-of-arrival distribution in the interacting case, 
both formally and numerically. This exam has yielded the result that, among these three, there is only one. Hi, in fact 
the simplest, that satisfies the correspondence principle in the sense of recovering the classical expression for the time 
of arrival when the effects of non-commutativity of the operators involved may be neglected. The other two proposals 
provide the correct classical limit in certain cases, but not in general. The numerical analysis of these distributions 
supports these formal considerations. 

In order to make an adequate comparison between these three proposals, it has proved convenient to write them in 
the unified formalism of the "crossing states". The formalism itself suggests further generalizations, by considering 
alternative crossing states. An open question is whether the crossing state formalism is indeed the most adequate 
one for the description of times of arrival, or whether other presentations are more suitable. What is assured is both 
that the crossing state formalism guarantees covariance and positivity, and its power for comparing widely diverging 
previous proposals, this being a property that we would expect to hold more generally. 
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